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Abstract 

In this paper, we study mixing and large decoherence in continuous-time 
quantum walks on one dimensional regular networks, which are constructed 
by connecting each node to its 21 nearest neighbors(/ on either side). In 
our investigation, the nodes of network are represented by a set of identical 
tunnel-coupled quantum dots in which decoherence is induced by continu- 
ous monitoring of each quantum dot with nearby point contact detector. To 
formulate the decoherent CTQWs, we use Gurvitz model and then calcu- 
late probability distribution and the bounds of instantaneous and average 
mixing times. We show that adding links to cycle network, in appearance 
of large decoherence, decreases the mixing times. Also, the mixing times 
are inversely proportional to squared parameter /, but are linearly propor- 
tional to the decoherence rate. Moreover, our analytical solutions justify 
the mentioned results for cycle in |18j and for long-range interacting cycle 
in [221. 
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1 Introduction 



The classical random walk is one of the most basic concepts of probability. 
When formulated with a finite number of states, it consists of randomly 
choosing paths between vertices on a graph, and many randomized algo- 
rithms for computing are such walks. Random walks have a clear appeal to 
both physicists and computer scientists: they are a powerful tool for both de- 
scribing physical phenomena and writing algorithms. Not surprisingly then, 
the quantum mechanical analog of random walk(quantum walk) can be de- 
veloped from both physical and computational directions as well. Quantum 
walks were formulated in studies involving the dynamics of quantum diffu- 
sion [Ij, but the analysis of quantum walks for use in quantum algorithms 
was first done by Farhi and Gutmann [2j. Quantum walks as noted earlier, 
are primarily of two types. Depending on the way the evolution operator 
is defined, they can be either discrete-time quantum walks(DTQWs) [3] or 
continuous-time quantum walks(CTQWs) [2|. In the CTQW, one can di- 
rectly define the walk on the position space, whereas in the DTQW, it is 
necessary to introduce a quantum coin operation to define the direction in 
which the particle has to move. Unlike the classical case the DTQW, because 
of the coin space, does not reduce to the CTQW when we let the time step 
between repeated applications of the unitary operator tend to zero. In recent 
years, the DTQW has been studied on random environments [4J, on quo- 
tient graphs [5j , in phase space [6j and for single and entangled particles f7] . 
Also, the CTQW has been considered on star graph [Hll], direct product of 
cayley graphs |lOj , quotient graphs [H] , odd graphs [12] , trees [HI E] and 
ultrametric spaces [15]. All of these cited articles have focused on the co- 
herent quantum walks, i.e. we have a closed system without any interaction 
with any environment. But in any realistic quantum system, e.g. quantum 
computer, interaction with the environment cannot be entirely eliminated. 
Investigations on the decoherent QW were given in [161 (171 dHJ (El [20] . Also, 
the effect of decoherence on the mixing time of CTQW has been provided 
on cycle [ISj, long-range interacting cycle \21\ I22j and one-dimension(lZ?) 
regular network [23]. Long-range interacting cycle is constructed by con- 
necting all the two nodes of distance m on the cycle graph [24] while ID 
regular network is a cycle graph in which every node is connected to its 21 
nearest neighbors(/ on either side) [25]. Figs. 1(a) and (b) show sketchs of 
LRIC with N = 8 and m = 3 and ID regular network with = 8 and 
/ = 3, respectively. It is clear that ID regular network can be regarded 
as the extensions of the long-range interacting cycle network in which the 
interacting parameters are m = l,2,-- - ,1. 
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(a) (b) 



Figure 1: Examples of LRIC with N = 8, m = 3(left) and ID regular 
network with N = 8 and I = 3(right). 

In [23], authors proved that the mixing time of CTQW on ID regular 
network, for the small rates of decoherence, is inversely proportional to rate 
of decoherence, but is independent of parameter 1(1 > 2) (due to losing of 
quantum interference). In this paper, we provide CTQWs on ID regular 
network in the appearance of large decoherence. We show that the bounds 
of instantaneous and average mixing times are proportional to decoherence 
parameter but, unlike the small decoherene case, are inversely proportional 
to the sum of squared interacting parameters(m). Moreover, our analytical 
results for I = 1 and for any m ^ 1 are in agreement with the results 
mentioned in [18] and [22], respectively. 

Our paper is structured as follows: In Sec. 2, we provide a brief summary 
of the main concepts and of the formulae concerning CTQWs and give the 
exact solutions to the transition probabilities on ID regular network. Sec. 
3 present the decoherent CTQWs on ID regular network. We assume that 
the decoherence rate F is large and calculate the probability distribution in 
Sec. 4. In Sec. 5, we obtain the lower and upper bounds of instantaneous 
and average mixing times. Conclusions and discussions are given in the last 
part. Sec. 6. 

2 CTQWs on ID regular network 

Every network can be considered as a graph whose nodes are connected in 
a similar manner. From an algebraic point of view, to each graph there is 
a corresponding adjacency matrix A = (Aij), which is a discrete version 
of the Laplace operator. The non-diagonal elements Aij equal 1 if nodes i 
and j are connected by a bond and otherwise. The connectivity of node i 
can be calculated as a sum of matrix elements Zi = Aij . The Laplacian 
operator is then defined as -L = Z — A, where Z is the diagonal matrix 
given by Zn. = ZiSik- It is worth underlining that, being symmetric and 
non-negative definite, L can generate both a probability conserving Markov 
process and unitary process. Thus, the Laplacian operator can work both as 
a classical transfer operator and as a tight-binding Hamiltonian of a quantum 
transport process [26 ] [27 ] [28]. 
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The continuous-time random walks(CTRWs) are described by the fohowing 
Master equation [29] : 

d ^ 
1=1 

being Pkj{t) the conditional probabihty that the walker is on node k at 
time t when it started from node j. If the walk is symmetric with a site- 
independent transmission rate 7, then the transfer matrix T is simply related 
to the Laplacian operator through T = —jL{m the following we set 7 = 1). 
The quantum-mechanical extension of the CTRW is called continuous-time 
quantum walk(CTQW). The CTQWs are obtained by identifying the Hamil- 
tonian of the system with the classical transfer matrix, H = —T [21 [271 130] ■ 
The states representing the walker localized at the node j, span the 
whole accessible Hilbert space and also provide an orthonormal basis set. In 
these basis the Schrodinger equation is 

i^^\k) = H\k), (2) 

where we set m = 1 and h = 1. The time evolution of state \j) starting at 
time is given by \j,t) = U{t)\j), where U{t) = exp[—iHt] is the quantum- 
mechanical time evolution operator. Therefore, the behaviour of the walker 
can be described by the transition amplitude Q.k,j{t) from state \ j) to state 
I k) , which is 

a.jit) = (fc|e-^*|i). (3) 
Prom Eq. (2), the a^jit) obeys the following Schrodinger equation: 

d ^ 

-rrakjit) = -i^Hkiaij{t). (4) 

1=1 

Note that the squared magnitude of transition amplitude provides the quan- 
tum mechanical transition probability ■Kkj{t) = \akj{t)\'^. 
To get the exact solution of Eqs. (1) and (4), all the eigenvalues and eigen- 
vectors of the transfer operator and Hamiltonian are required. We denote 
the nth eigenvalue and eigenvector of A by En and |g„), respectively. Now 
the classical probability is given by 

N 
n=l 
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and the quantum probability can write as 

N 

^kAt) = E ^~''^''"~'''Hk\Qn){qn\j){k\qi){qi\j), (6) 

n,l=l 

In the following, we focus on ID regular networks and study CTQWs on 
them. ID regular network is constructed by connecting every node to its 2m 
nearest neighbors(?Ti on either side), see Fig. 1(b). This network provide a 
good model to study various coupled dynamical systems, including biological 
oscillators [31) . Josephson junction arrays [32], synchronization [33j, small- 
world networks |34j and many other self-organizing systems. The Hamilto- 
nian of the system is given by [23] 

[ -2/, ifi = j; 
Hij = < 1, if i = j ±m,m e [1,1]; (7) 
[ 0, Otherwise. 

This Hamiltonian acting on the state \ j) can be written as 

I 

H\j) = -{2l + l)\j)+ J2 + (8) 

m=—l 

which is the discrete version of the Hamiltonian for a free particle moving 
on a lattice. It is well known in solid state physics that the solutions of the 
Schrodinger equation for a particle moving freely in a regular potential are 
Bloch functions [35l |36]. We denote the Bloch states by \^n) and then the 
time independent Schrodinger equation can be written as 

H\^n)=En\^n). (9) 

The Bloch state \^n) can be expressed as a linear combination of the states 
|j) localized at nodes j, 

N-l 

^ i=o 

Because of periodic boundary conditions, we have <I>„(A^) = <I>„,(0), where 
^n{x) = This restricts the 0„-values to 9n = ^jp^, where n = 

0, 1, . . . , — 1. We set Eq. (8) and (10) into (9) and obtain the eigenvalues 
of system as [23] 

I 

En = -2l + 2Y,cos{j9n). (11) 
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Figure 2: Fig. 2(a) shows point contact detector j monitoring the electron 
in dot j + 1 and Fig. 2(b) shows point contact detector j when electron is 
placed in dot j. Source and drain reservoirs are kept at zero temperature 
with chemical potentials Fij and Frj, respectively. Ej is the on-site node 
energy. Qirj and Cli^j are transmission probability of the detector for empty 
dot and for the occupied dot. 

Thus Eq. (5) and (6) can be written as 

n 

,r,^^.(t) = l^^e'tiE„-EO^-iik-mn-i)'^_ (13) 

n,l 

3 The Decoherent CTQWs on ID regular network 

In this section, we take into account the effects of decoherence in quantum 
walks on ID regular network. For this aim, we use the realistic physical 
model of the network in which decoherence induced by continuous moni- 
toring of each network node with nearby point contact(PC) detector. In 
this model, our assumptions are as follows: each node is represented by a 
quantum dot which is continuously monitored by an individual point con- 
tact(PC), the walks are performed by an electron initially placed in one 
of the quantum dots, identical PCs are placed far enough from QDs so 
that the tunneling between them is negligible. Coulomb interaction between 
electrons in QD and PC is taken into account and all electrons are spin-less 
fermions |37] . Such a set up is shown schematically in Fig. 2, where detector 
j is represented by a barrier, connected with two reservoirs at the potentials 
Fij and Frj. The transmission probability of the barrier varies from ^irj 



6 



to ^'irji depending on whether or not the quantum dot is occupied by an 
electron. 

In the following, we want to write Hamiltonian for the entire system. We 
consider simple continuous-time quantum walks are defined over an undi- 
rected graph with nodes in which each node is labelled by an integer 
i G [0, — 1]. These walks can be well described by the following Hamilto- 
nian [MlEHj: 

(14) 

= J]A.,(t)(|i)(j| + |j)(i|)+^i^,-(t)|j)(i|. 

ij j 

Here, each node j corresponds to the quantum state |j) = cj|0) in which 
\j) denotes the state where the particle is located at node j. The two 
terms correspond to a hopping term with amplitudes Ajj(t) between nodes, 
and on-site node energies Ej{t), both of which can depend on time. For 
the sake of simplicity, we drop all the on-site energies and assume that 
hopping amplitudes(Ajj(t)) between connected sites to be constant. Also, 
we renormalize the time, so that it becomes dimensionless [37]. Thus, Eq. 
(14) reduces to 

N-l I 
j=0 z=l 

Now, we study the electron transport process in the point contact j. The 
point contact is considered as a barrier, separated two reservoirs (the source 
and drain) . All the levels in the source and drain are initially filled up to the 
Fermi energies. We call it as the vacuum state |0). Thus, the Hamiltonian 
of j-th point contact can be written as 

Hpc,j = Y + Y ^^j^lj^^j + Y ^irjialjcirj + aljttij), (16) 

I r l^r 

where a} j{0'ij) and j(ar,j) are the creation(annihilation) operators in the 
left and right reservoirs respectively, and ^irj is the hopping amplitude 
between the states Eij and Erj in the right and left reservoirs. 
The interaction between the detector and the measured system is described 
by Hint- The presence of an electron in the left dot results in an effective 
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increase of the point contact barrier (fi^rj = ^ir,j + ^^irj)- Therefore, the 
interaction Hamiltonian can be written as 

Hint = ^ ^^IrjCjCjialjClrJ + al jCLij). (17) 
l,r 

For simphcity, we assume that the hoping amplitudes are weakly dependent 
on states Eij and Erj, so that Jl/rj = ^^irj = and Fij{Frj) = 
Fi{Fr). Gurvitz in [ID] showed the appearance of decoherence leads to the 
collapse of the density matrix into the statistical mixture in the course of 
the measurement processes, thus the evolution of the reduced density matrix 
traced over all states of source and drain electrons is given by Bloch-type rate 
equations. The time dependent non- unitary evolution of reduced density 
matrix p{t) = \^{t)){^{t)\ in the Gurvitz model is given by [23] 

/ 

■SiPjM't) = -ji X] {Pj+m,k - Pj,k+m)] - r(l - 6j^k)Pj,k 
m=—l 

I 

— {Pj+m,k — Pj,k+m + Pj-m,k — Pj,k-m)] — r(l — 6j^k)Pj,k- 

4 Large Decoherence 

In [23], authors studied the effect of small decoherence (TiV ^ 1) in CTQWs 
on \D regular networks. They showed that the instantaneous mixing time 
upper bound and the average time mixing lower bound are independent of 
parameter /(/ > 2), but are proportional to the inverse of decoherence rate. 
Here, we want to study CTQWs on these networks in appearance of the 
large rates of decoherence. For this purpose, we use Gurvitz model and 
focus on the elements of matrix p{t). Based on the initial conditions, the 
non-zero elements appear only along the major diagonal. Firstly, we rewrite 
Eq. (18) for the elements of major diagonal and also minor diagonals whose 
distances of the major diagonal are lesser than /. By dropping terms that 
are smaller than l/F, we have 
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Pj,j^^ ~ iV/, iPj+m,j Pj,j+m ''r Pj-m,j Pj,j-m)]j 
m=l 

^ p'jj+i{t) = -i[pj+i,j+i - Pj,j] - Tpj-j+i, 

P'j,j+2{'t) = -\[Pj+2,j+2 - Pj,j] - rpjj+2, 
^ P'j,j+lit) = -{[pj+lj+l - Pjj] - Tpjj+i, 

For the simplicity sake, we use the following definitions 

~ Pj,j 
= Pjd+i " Pj+i,j^ 
fj = PjJ+'^ ~ Pj+'^,j' 

1i = PjJ+i - Pj+iJ- 
Then the above difference equation system can be written as 







-1 - fj + /7-2 + (I3+ Qj-l 




-|[aj+i - aj 


-Tdj, 


< 


- 0-3. 











Differentiation of the above equation gives 
d'; = -m^,-a',]-Tdf,, 

/;=-i[«;+2-a;]-r/i, 

We can guess the following solutions for the above equations 



AT-l 



E 27rijfc 



-7fc* 



N-1 



k=0 

N-1 



J- \ ^ rn -'tut 



k=0 

N-1 



(23) 



E 27Tijk 



-ikt 



k=0 



k=0 



in which jj^, Aj^, D^,--- and are the unknown quantities. We set 
these solutions into Eqs. (22) and get 



7feAfc + |[L>fc(l - eT) + Ffe(l - eT) + • • • + Q;i^(l - e^^)] = 



Akiii-e^ + 1)] + Dki-fk - r) = 0, 



4:7vik 



(24) 



+ 1)] + Fki^k - r) = 0, 



I Afe[i(-e^ + 1)] + Qkhk - r) = 0, 

It is evident that there are nontrivial solutions when the determinant of the 
coefficients matrix is zero. 

(7. - r)('-^)[7.(7. - r) + ^(-n^(^) + sin^(^) + " " " + -n^(^))] = 0(25) 
Therefore, there are four values for 7*; as 



7fc,o = 0, 

7fc,i = r, 



Ik 



7fc,2 
7fc,3 



2T ^ 



m=l 



N 



(26) 



1 \ - . n.Trfcm, 
— > sm'^( ) 



m=l 



The general solutions of Eqs. (22) are 
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N-l 



k=0 

N-l 

(27) 

Af-l 
k=0 



N-l 
k=0 

By the initial condition aj{0) = 6jfi and dj = 0, for j = 0, ••• ,iV — 1, 
and replacing the constant coefficients into the others, we get 

Ak,i ^ 0, Ak,2 - ^™^(-^)' ^k,3 = 1, 

m=l 

L>fe,i~0, L»;,,2~^sin(^)exp(^), Z^^^3~=^sin(^)exp(^), ^^g) 
~ 0, Fk,2 ^ f sin(2^) exp(2f^), ^^,3 ^ ^ sin(2^) exp(2^), 

I ^ 0, Qk,2 ^ ^ sin(^) exp(^), Qk,^ ^ ^ sin(^) exp(^). 

Note that the probability distribution P{t) of the quantum walk is specified 
by the diagonal elements of p{t), that is Pj{t) = aj{t). For large rates of 
decoherence, this distribution reduces to 

fc=0 

Since CTQWs are symmetric under time-inversion, the above distribution 
does not converge to any constant value. 
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5 Mixing time 



The rate of convergence to probability distribution of CTQWs can expressed 
in terms of many quantities, but the one used commonly is mixing time. 
There are two different definitions of mixing time in the literature: the in- 
stantaneous mixing time and the average mixing time. To define the notion 
of the mixing time, we use the total variation distance between distributions 
di and d2 as \\di - d2\\tv = Ej " d2{j)\- 

(a) Instantaneous mixing time 

The instantaneous mixing time is defined as the first time instant at which 
the probability distribution of the walker position is sufficiently close to 
uniform distribution [38], i.e. 

7V-1 ^ 

Tins = min{T : ^ \\Pj(t) - — \\tv < e}, (30) 
j=0 

Adding a summation over j to Eq. (29) gives us the total variation 
distance required to obtain the mixing time: 

Af-l N-1 N-1 I J 

E Mt) 4i = E 1^ E -p[g E - ^1 

j=0 j=0 k=0 m=l 

j=0 k=l m=l 

Lower bound 

To find the lower bound, we apply the following way: 

Firstly, we use the term j = and then use the terms k = 1,N — 1 

j=0 k=l m=l 

2 ~2r E (^) 




m=l 
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Based on the instantaneous mixing time definition (Eq. (30)): 



-e rn=l =e. (33) 



Thus, the lower bound of mixing time has the form 



For iV > 1 



2T 2 

^ins,lower — j llli^ ). (34) 

m=l 



2rA^ 2 



2 2 

m 



TO=1 

Note that for / = 1, we have 



_ 2riV2 .2 

J-ins,lower — n ^■'^VlTj^/ l^Oj 

vr^ iVe 



which is coincide with [18] 's result for cycle. 

Also, it is evident that for any m 7^ 1, relation(36) gives us the mentioned 
result in [22] for LRICs. 



Upper bound 

To obtain the upper bound of instantaneous mixing time, we use the follow- 
ing way: 

E l-.W - ^1 = E I^E -■'ig E sin^(t^)lcos(M5)| 

i=o i=0 k=l m=l 

N-lN-l , I 



4 E E -Pig E 

j=0 fc=l m=l 
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Because of sin(x) > ^ for < x < | ^4T], we have 
^-1 1 2 f"^/'' -t ' nkm 

j=0 i=0 fc=l m=l 



< 



< 



< 



j=0 k=l m=l 



j=0 fc=l m=l 



2 



2^exp[— 

j=0 A:=l m=l 



that in the third inequahty, we used of the relation k"^ >k for k>l. 
some algebra, we obtain 



N-l 



j=0 ^ A^2m^ 



m=l 

According to the instantaneous mixing time definition(Eq.(30)): 

2 

; = e. 



"^P[ r ^ ^1 " ^ 

m=l 



Therefore, the upper bound of instantaneous mixing time is 

TA^^ /2 + e 

-'■ ins , upper — ^ /• 
m=l 

For / = 1 (cycle), we obtain 

_rw2 2 + e 

-I- ins , upper — „ ) 
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and for any m 7^ 1 (long-range interacting cycle), we have 

_ riV2 2 + e 
2(1 + m^j e 

which are the same results mentioned in |181 [22] . 
(b) Average mixing time 

The average mixing time is based on the time-averaged probability distri- 
bution, that is 



1 



T 



P{j,T) = - P{j,t)dt. (44) 

We can define the average mixing time as the number of time steps required 
for the average distribution to be e-close to the limiting distribution [42) . 
i.e. 

Tave = min{T : \\P{j,t) - < e}. (45) 

Firstly, we want to calculate P{j, T) by Eq. (29) 

■'^ k=0 m=l 

1 [ ^ ■ 2/^^"iM jk (46) 



m=l 



TN ^ I 1 

m=l 

The total variation distance between the uniform distribution and the time- 
average distribution of the decoherent quantum walk is given by 
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N-1 _ ^ N-l 2p ^-1 



krr 

w 



N-1 ^ ^1 2r 2r ' ■ 1 



sin (- 

m=l 



N-1 7V-ll-«^P[-^E«i^'(^^l'^'' 



'TN ^ I 

sm (- 

m=l 



j=0 k=l V r.in^r^^"^ 



Lower bound 

The lower bound of average mixing time can be derived as follows: 



TV ' - TiV ^ i 

E--'(- 

m=l 



> sm ( 

N 



(48) 



1 — exp — — > sm 1—:^) 

^ m=l 



TN I 

E (^) 

m=l 

where the first inequality uses the term j = only and the second inequality 

uses the terms fc = 1, A'' — 1. 

Assuming T/T ^ 1 (which is consistent with our other assumptions, N ^ \ 
and FA^ 1, it requires T ^ N) gives 



Af-l 

\P(i.T'\ - 

N 

TN siu^r 

N 



1 4r 

Y.\Pi^^T)--\> (49) 



TAr^sin2(_) 



m=l 
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(50) 



Therefore, the lower bound of average mixing time is 

_ 4TN 

J-ave,lower — „ "Z- 
Em=l " 

Upper bound 

Now, the upper bound of average mixing time can be obtained as follows 

N-l , K-l „„ K-1 1 - '=''P[-^ 2^ ™ ("i^'l'^ 

=i:i|^E r^^, 1 

m=l 

1 — exp — — > sm {— -— ) cos( — -— ) 
^ 2r N N 

m=l 



N-l 2p ^-1 



-y\—y- 

m=l 

r T J-^ . r, Trkm. 
^_^l_exp[-— ^sm (^ 

m=l 



N-l AT-il-exPh^ 2^sm 

m=l 



iv 

-) 



< 



[iv/2]l-exp[-2p 2^sm (— 

\ ^ m=l 



TN ^ I , 



■)] 



m=l 



2T ^ A:2m2^ 
[7V/2]l-cxp[-— 

<iy "'=^ ' 

fe=i fc 

m=l 
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We assume that F/T » 1 and obtain 



m=l 

< —, C(2) < -- 



r7i=l m=l 



where C is Riemann zeta function [43j. 

Thus, the upper bound of average mixing time is 



m=l 

Therefore 



ln(— ) < 

m=l 



ATN riV2^2 



2 2 "m^^.a^'e ^ r 2 



(52) 



(54) 



Comparing of the both mixing times shows that the instantaneous mixing 
happen earher than the time-average mixing. The above relations show that 
the bounds of instantaneous and average mixing times in the decoherent 
CTQWs on cycle decreases with adding the newly edges. Moreover, the 
bounds of mixing time have linear dependence on decoherence rate, while one 
observes inverse dependence on the sum of squared interacting parameters. 



6 Conclusion 

We considered the continuous-time quantum walks on one-dimension regu- 
lar network under large decoherence T 3> 1. For this, we used an analytical 
model developed by Gurvitz [40] and calculated the probability distribution 
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analytically. Then we calculated the lower and upper bounds of instanta- 
neous and average mixing times as 



Thus instantaneous mixing time is shorter than average one. The both mix- 
ing times are linearly proportional to the decoherence rate and are inversely 
proportional to the sum of squared interacting parameters (m). We found 
that adding shortcuts to cycle network decreases the mixing times. More- 
over, our analytical results for I = 1 and any m ^ 1 are agreement with the 
mentioned results in [18} I22j . 
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